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1. Introduction 

The proposals of A. B. Zamolodchikov for the scaling region of the Ising model in a 
magnetic field 0] have stimulated a renewal of interest in two-dimensional integrable field 
theory. To think of a model as a perturbation of a conformal field theory by certain of its 
relevant operators allows much non-perturbative information to be extracted, starting with 
the existence of the conserved charges which guarantee integrability, and often leading on 
from there to hypotheses for an exact S-matrix and mass spectrum, hypotheses which in 
many cases have now withstood independent checks. 

One such check is to 'complete the circle', reconstructing information on the ultraviolet 
limit from the infra-red data contained in the hypothesised exact S-matrix. If this matches 
with the content of the initial conformal field theory, a very non-trivial confirmation of 
the whole approach has been achieved. Better still, non-perturbative information on the 
whole renormalisation group trajectory between ultraviolet and infra-red can often be 
extracted. The machinery employed is known as the Thermodynamic Bethe Ansatz (TBA) 
0, first used in this (relativistic) context by Al. B. Zamolodchikov ||; subsequently these 
techniques have been developed by many others. 

The method is straightforward if the scattering is purely elastic 0-0 , that is if the 
S-matrix is diagonal in a suitable basis. The thermodynamics of the system in a box can 
be reconstructed from a set of coupled non-linear integral equations (the TBA equations), 
and these equations follow immediately from the logarithmic derivative of the S-matrix 
itself. If the box is small enough, UV effects dominate and (to mention one example) the 
central charge of the original conformal field theory can be recovered. Interesting features 
have been observed in the study of some of these systems of equations, pointing to the 
possibility of further exact results 0: the TBA equations for models described by the so- 
called minimal ADE affine Toda S-matrices (see below) can be written in a universal way, 
depending only on corresponding non-affine Dynkin diagrams. 

Unfortunately, the majority of cases involve non-diagonal S-matrices. To derive TBA 
equations directly from such scattering data, the so-called 'colour transfer matrix' should 
be diagonalised |§ , a forbidding problem - requiring the higher-level Bethe ansatz (see, 
for example, ref. || ) — in many cases. Furthermore, many perturbations of conformal field 
theories lead to flows, not to massive scattering in the IR limit, but to a non-trivial fixed 
point so that the RG flow interpolates between two conformal field theories. The ^13 flows 
between minimal models |1CJ] provide the first examples of this phenomenon. Despite some 
suggestive proposals | ]TT]JT^ ] the S-matrix approach is much more delicate in such cases, 



owing to the well-known difficulties of treating massless particles in 1+1 dimensions. 
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An alternative strategy, developed by Al. B. Zamolodchikov in a series of papers 



f8|, |TT| , |T3|| , is to start with an educated guess for the relevant system of TBA equations, 
and then to check that they correctly reproduce a reasonable number of parameters typ- 
ical of the UV theory, such as the central charge and the dimension of the perturbing 
operator. Once it has been established to general satisfaction that the proposed TBA is 
indeed associated with the perturbed conformal field theory under discussion, the RG flow 
of the vacuum energy (and possibly also of the energies of some of the excited states |T^] ) 
can be followed numerically, or IR data can be reconstructed analytically, allowing the IR 
limit to be identified - perhaps with some other conformal field theory. 

In this way Al. B. Zamolodchikov was able to propose TBA systems for the flows of 
minimal models perturbed by </>i3, the least relevant operator. Rather remarkably, these 
TBA systems can also be encoded on graphs (in fact the A n Dynkin diagrams). In the 
diagonal TBA, energy terms in the equations, one for each particle type, were naturally 
attached to the nodes of the corresponding Dynkin diagram, as will be explained further 
below. Here, the structure is a little different, with all but one or two of the energy terms 
being set to zero. The correct distribution to take depends on the sign of the perturbing 
parameter, and has a dramatic effect on the predicted IR behaviour: for one choice the 
flow is to a massive theory of kinks || , while the other gives the flow to the next- lowest 
minimal model |llj- Subsequent work|0|], further modifying the distribution of energy 
terms, gave ansatze (with identical diagrammatic structure) for TBA systems describing 
the flows between su{2) k x su(2) 1 j su(2) h+l coset models. 

Proceeding further, it was natural to look for TBA systems to describe analogous 
perturbations of g( k >xg">/g( hJrl > coset models, for which </>i3 becomes (j)(id,,id,adj) ■ These 



were given for g = A n , I = 1 in [15|, and then for the general (simply-laced) case in [16] . 
In this latter reference a universal form of the general TBA equations was given, allowing 
the generalisation of the diagrammatic structure found in |7] to be identified: the relevant 
graph is a 'product' of the Dynkin diagram of g (which had already been seen to be relevant 
for the diagonal Toda-type S-matrices) with an Ak Dynkin diagram. Otherwise stated, 
the graph is formed by taking a tower of k copies of the g Dynkin diagram, with links in 
the vertical direction joining the replicated nodes. Presumably this replication encodes the 
colour structure of the non-diagonal S-matrices for these models RITA , although the details 
remain to be worked out. 

For the minimal models (g = su(2)), there is another intriguing result. Al. B. Zamolod- 
chikov recently proposed |1£| a simple diagonal scattering theory (the so-called 'staircase 



model'), from which a TBA system followed in the standard way. However, numerical 
study of the solutions of this system revealed a structure much richer than any previously 



observed, with hoppings of the vacuum energy (which is proportional to the effective cen- 
tral charge) strongly suggestive of an RG flow for an underlying field theory that passes 
by a sequence of c < 1 (minimal model) fixed points. 

In fact, Zamolodchikov's S-matrix is an analytic continuation of that for the A\ affine 
Toda field theory. This suggests a natural generalisation to the other affine Toda field 
theories, and our aim in this paper is to present arguments supporting the claim that such 
S-matrices yield new staircase models which hop between the g( k >xg( 1 >/g( k+1 ) (W g ) coset 
conformal theories. The above-mentioned structure of a tower of Dynkin diagrams for the 
TBA systems found in [Tnj emerges as a natural consequence of the block structure of the 



affine Toda S-matrices. 

While we were completing the writing of this paper, independent work by Martins 
appeared |H| , proposing the same generalisation of |TB| as we give here, and backing this 
up with a detailed numerical study of the solutions of the TBA system. This seems to 
us to complement our rather more analytic approach, and we therefore refer the reader to 
|T9t1 for numerical confirmations of the various approximations and hypotheses to be made 
below. 

The paper is organised as follows: in section 2 we recall the necessary features of 
the affine Toda S-matrices, and their relationship with the root systems of the simply- 
laced Lie algebras. The connection between the 'minimal' and 'full' affine Toda S-matrices 
is described, and this leads to a key relationship between their logarithmic derivatives. 
Section 3 presents the relevant TBA equations, and discusses the 'locality' structure of 
these equations as a function of the coupling constant. For a suitable analytic continuation 
of the coupling constant, the system becomes highly non-local and involves two scales 
(rather than the more usual one) , the interaction of which is responsible for the intricate 
hopping structure of the solutions. Effective central charges are predicted for various 
asymptotic limits of these two scales, providing evidence that the RG behaviour being 
described is indeed a roaming trajectory visiting a sequence of M^-minimal models. At 
the risk of repeating material that may be well-known to experts in the field, we have 
decided to make this part of the discussion reasonably explicit. Finally, we give a few 
speculations as to the physical meaning of the structures observed. 

2. The affine Toda S-matrices 

For real values of the coupling constant, the simply-laced affine Toda theories seem 
to have a number of particularly simple properties p0|-|23[| . There are r different types of 
particles, r being the rank of the associated Lie algebra. Particles of different types are 



distinguished one from another by conserved charges of non-zero spin, constraining the 
S-matrix to be diagonal — the scattering theory is purely elastic. Furthermore (and in 
contrast to the situation for non simply-laced algebras) the mass ratios remain fixed at 
their classical values, and do not pick up any coupling-constant dependence via quantum 
effects. As a result, the full S-matrix elements factorise into a product of two pieces: a 
minimal factor, independent of the coupling constant, which contains the physical-strip 
poles and thus all information on mass ratios and bound-state structure; and a second 
'Z-factor' incorporating the coupling-constant dependence. This piece simply introduces 
some zeroes into the physical strip, ensuring that as the coupling constant tends to zero 
the S-matrix tends to the identity, but changes neither the bound-state structure nor the 
resulting bootstrap equations. The situation for non simply-laced theories is markedly 
more complicated [[23] , and will not be discussed below. 

Despite their apparent independence, both factors must satisfy the same set of boot- 
strap equations, equations which follow from the pole structure of the minimal piece. As 
a result, their functional forms are very similar, a fact which will be relevant below. This 
is most clearly seen when the S-matrix elements are written as a product of the functional 
building blocks {x}b employed in ref. P2|j . First, define elementary blocks (x) by 

(?) (0) = )l 2h [ , (2.1) 

sinh(f-^) 

where h is the Coxeter number of the algebra involved. (Other authors use the notation 
f a (d) for (ha) (9); including the factor of h in the notation, as here, has the advantage that 
the parameter x then turns out to take integer values.) The bigger block, incorporating 
the coupling-constant dependence via the function B((3) = 2/3 2 /((3 2 +4tt), is 

r-, (X-1)(X + 1) 

{X}B ~ ( X -1 + B)(x + 1-B) ■ (2 - 2) 

The factorisation of this block into minimal and coupling-constant dependent pieces is 
clear, as is the relationship between the two factors. To express this precisely, define a 



symmetrised shift operator pljl , acting on functions /(#), by 

(TZ y f)(9) = f(9-my/h)f(9+my/h). (2.3) 

Then 

K y (x) = (x-y)(x + y), (2.4) 

and so 

{ x j b = 15 — rr- 2 - 5 ) 

H\-B \x) 
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Thus the only difference between the two factors in the building block is a reciprocation 
and a change in the shift. This property is inherited by the full S-matrix elements, though 
it can sometimes be hidden due to the relations (0) = — (h) = 1. 

From the point of view of the bootstrap equations, it is perfectly consistent to consider 
the minimal part on its own, as removing the Z-factor leaves the physical poles and the 
signs of their residues unchanged. Such S-matrices are products of the corresponding 
minimal blocks, 

{x} min = (x - 1) (x + 1) = n x ix) . (2.6) 

However, the physical interpretation is very different: TBA calculations [^-^J have es- 
tablished beyond reasonable doubt that these S-matrices correspond to perturbations of 
g^'xg^'/gC 2 ' coset conformal field theories. 

Although the details are not essential for the main development, it is worth recalling 



how these S-matrix elements can be written in a uniform way for all the ADE series |25|] . 
The construction uses properties of a Coxeter element w taken from the Weyl group of the 
relevant root system $, each particle species a being associated with an orbit T a C $ of 
w. Letting r+ = T a n $ + denote the intersection of such an orbit with the positive roots, 

ST = II W», 0) + l}^ ■ (2.76) 

£er+ 

Here, A a is the fundamental weight associated with the a th spot on the Dynkin diagram, 
4> a =(l— w~ 1 )X a is a particular root in the orbit T a , and 27ru(0 a , (3)/h is the angle between 
the projections of (p a and (3 into the e 27 ™' ^-eigenplane of w. For more explanation see [55 
but really it is only the general form of ( |2.7[ ) that will be relevant to what follows. 



At the level of the S-matrices, the factorisation seen for the blocks in ( |2.5| ) can now 
be made explicit by introducing 



With this definition 
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= K-iS ab 
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(2.9) 

Note, the sign in (|2.8|) is left undetermined by the requirement that ( |2.9|) should hold; 
indeed S^ b can be thought of as a skewed 'square root' of S™ b in . As long as the function 



is used merely as a book-keeping device, this ambiguity will not cause any problems. 
However, it is tempting to speculate about a physical interpretation, in which case the 
sign might perhaps be important and a further subtlety should be considered: namely, 
whatever the sign choices made in (|2.8), the functions S^ b will in general fail, by a sign, 



to satisfy some of the bootstrap equations holding for S^ n and S^iP- The origin of these 
signs can be seen by re-examining the proof, given in ||25|| , that the bootstrap equations 



hold for the functions defined in ( |2.7|) . Since the bootstrap equations often involve shifts 
in arguments more general than the symmetrical operation defined in ( |2.3| ), the first step 
must be to write (|2.7| ) using functions which transform into themselves under the single 
shift T y , defined by 

(T v f)(0) = f{0 + i*y/h). (2.10) 

This property does not hold for (x) defined in (pTT|), but decomposing it as (x) = 

( x ) + / ( — x ) + 5 where 

(x) + (6) = sinh^ + S , (2.11) 



results in a building block which satisfies 

T y (x) + = (x + y) + . (2.12) 

Similarly {x}_|_ can be defined by replacing (.) by (.) , in either (|2.2| ) or ( |2.6|) , in terms of 
which ( |2.7|) becomes 

S ab = [IWi.^ + l}^' (2.13) 

Per b 

(The root /3 now runs round the whole orbit Ti,, rather than just the positive half as in ( |2.7| ) ; 
to see the equivalence of (|2.13|) and ( |2.7|) needs a few root system identities, which can be 



found in ||25| .) Now an intermediate step in the proof of the bootstrap equations involves 
re-expressing the product appearing in ( |2.13| ), and for this, the relation {x+2h} + = {x} + 
must be used. It is this property which fails by a sign when the same trick is attempted for 
SE, since (x+2h) + = — (x) + . Thus if the manipulations involve shifting an odd number 
of blocks (x) , , an extra minus sign will appear, and even the freedom of signs allowed 
by ( |2.8| ) is not enough to eliminate all these signs from the bootstrap equations (this can 
already be seen in the A^ case). There may be some interesting features hidden in these 
signs, but the question will not be pursued here. The signs in any event disappear once 
the logarithmic derivative has been taken, and so do not affect the TBA equations. 



this sign was missed in [16]. 



To discuss these logarithmic derivatives, set </>/(#) = — i-4log/(0). Following from 
this definition, 

tfg = 4>f + <^> ^f- 1 = -<t>f> <t>K y f = (Ty + T-y)4>f- (2-14) 

With the shorthands cj)^ b , <^ b in and <j)^ b defined in the obvious way, combining ( |2.9|) with 
( glD yields 

€ h = €T - 71-b^ - T B -x€ h - (2-15) 



For consistency with [IB], we write <p^ b as ?/> a &. Then, putting 

6 = -(l-B), (2.16) 

equation ( p. 15| ) can be rewritten as 

4&W) = <Pa b [n (o) - <M* + e ) - i> ab (e - e ). (2.17) 

This relation will be useful in the treatment of the TBA equations. 



3. The TBA equations 

The TBA equations for the affine Toda models have the standard form for a purely 
elastic scattering theory ||. Writing the mass ratios as fh a =m a /rni, a system living on a 
circle of finite length R depends on the scale through the quantity ^miR = e x , and the 
ultraviolet limit corresponds to sending x — > — oo. To avoid the complications caused by 
the fact that x will almost always be negative, this parameter will be swapped for y = —x. 
The TBA equations for the pseudoenergies e a (0), encoding all of the finite-size effects, are 



e a (9) = 2m a e-y cosh(#) - -L V ^ b * L b (9), (3.1) 

Z7T * — ' 
6=1 

where L a {6) = log(l + e~ 6a ^) and * denotes the rapidity convolution: 

/oo 
d6'(P(6-6')L(6'). (3.2) 

-oo 

The energy terms mentioned earlier are the functions 2fh a e~ y cosh(#); they can be attached 
to the nodes of the Dynkin diagram via the association P, |22|J2^ 1 between the Toda particle 



masses and the Perron-Frobenius eigenvector of the Cartan matrix. The effective central 
charge, defined to be —tt/6R times the finite size ground state energy E(R), is given by 

c(y, ) = —J2 d6m a e- y cosh(9)L a (9), (3.3) 

the #o-dependence of this function entering via the S-dependence of the kernel function 

To gain a general understanding of the solutions to Q3.1|) , the form of this kernel should 
be discussed. Since 

,„s 1 , (9 inx\ 1 , (9 %-kx\ ,„ tS 

wW = -coth I- + — -jooth (--— , (3.4) 



it follows from the definitions ( 2.76Q and Q2.8| ) (and the reality of the numbers u((fi a: f3) 



that for real 9 the functions <fi™ n (9) and ifi a b(6) are strongly peaked about 9 = 0, and 
decay like e - ' 61 ' outside an interval of size 0(1) centred at the origin. Thus so long as L(9) 
is reasonably well-behaved on the real #-axis, the values of (p min * L(6) and ift * L(6) are 
governed by the values taken by L(6') with 9' close to 9, and the convolution ( |3.2|) has 
a 'local' character. (To gain an intuition for this, it is helpful to think of replacing the 
functions (/> mm and ip by two delta-functions, each normalised to preserve the corresponding 
overall integral; indeed we shall be working in asymptotic regions where this seems to be 
a valid approximation to make.) For real values of the coupling constant (3, the shift 
#o in ( |2.17| ) is purely imaginary, and the above line of reasoning goes through for <p B as 
well, with the conclusion that the full TBA equation also has this local character. But 
to incorporate Zamolodchikov's roaming idea, and in particular to introduce another scale 
into the problem, it is necessary to be a little more general. This should involve relaxing the 
requirement that f3 be real, but on general grounds it might be expected that the S-matrix 
elements should continue to be real analytic - that is, S a b(9) should still be real for 9 on the 
imaginary axis. From the definition (|2.2|) of the basic block, any real value for B satisfies 
this condition, corresponding to (5 real if B is between and 2, and (3 purely imaginary 
otherwise. A little less obvious, but clear enough once ( |2.5| ) is considered, is that values of 
B such that (1— B)* = —(1—B) are also permitted (here, * denotes complex conjugation). 
That is, 1—B should be imaginary, and 9q real. In turn, this requires that /3/(2 v /7r) should 
lie on the unit circle in the complex plane. An appealing feature of the S-matrix elements 
for such values of (3 is that (unlike the situation for many imaginary values of (3) no extra 
poles are introduced into the physical strip, and so the interpretations of all the physical 
poles, bound states and bootstrap equations remain unchanged. It also connects with an 



interesting physical interpretation of the roaming behaviour, a speculative point that will 
be described in the concluding section. For now, the important point is that this analytic 
continuation changes the 'locality structure' of the TB A equations dramatically, as can be 
seen by using (|2.17| ) to rewrite ( |3.1| ): 

e a {9) = 2m a e- y cosh(#) 

'Til W * Lb{9) ' ^ ab * ( L ^- o) + L b (0+0o))] ■ (3.5) 

The mm term continues to 'see' Lb in the region 9' ~ 9, but for the ip term it is now the 
values taken by Lf, at 9 ± 9q, rather than at 9 itself, that are important. Furthermore, 
the equation now contains two scales, namely y and 9q. It is the interaction of these 
two scales which produces the staircase behaviour, the controlling parameter being their 
ratio, y/9o. If this ratio is kept fixed while y and 9q tend to infinity, then generally the 
effective central charge c(y, 9q) will tend to the central charge of one of the Wg-minimal 
models - precisely which one depends on the chosen value of y/9o. Different asymptotic 
behaviours are separated by integer or half-integer values for this ratio. If on the other 
hand y is increased while 9 is kept fixed, to study the behaviour towards the ultraviolet of 
one particular theory at one particular value of the coupling constant, then the changing 
value of y/9o results in a series of different values for the effective central charge, with 
a cross-over at each integer and half-integer. The rest of this section is devoted to some 
analytical arguments in support of these claims, leading to the expectation that the set 
of solutions to ( |3.5|) at given (sufficiently large) 9$ do indeed form a staircase, with the 
central charge of a g ( - k 'xg ( - 1 >/g ( - k+1 > coset conformal field theory at the k th step. 

As a first step, it seems reasonable to suppose that for fixed y the function L&(6>) is 



bounded. As soon as \9\ is large enough for the first (energy) term in (|3.5| ) to contribute, 
the exponential growth of the hyperbolic cosine will quickly swamp any influence of the 
convolution terms. Thus, for 9 ^> y or 9 <C —y, e a (9) is closely approximated by the 
functions £+ (9) or e~ (9) respectively, where 

et(9)=m a e-y ±e . (3.6) 

Correspondingly, L a {9) suffers a double-exponential decay beyond \9\ > y and is soon 
approximately zero. This information is already enough to see why the k=l fixed point 
is expected to control the effective central charge for y <C 9o/2. The reason is simply 
that for such a value of y, the V> a b term never contributes significantly to ( |3.5| ) - - for 
\9\ 3> y, the energy term dominates and all the convolution terms can be ignored, while for 

9 



\9\ < U, both \6 — 9q\ and \9 + 9q\ are larger than y, and hence Lb(9—9a) and Lj,(0+0o) are 
essentially zero, and the ifj a b term in ( |3.5| ) can again be dropped. The resulting 'reduced' 
set of equations is just that for S^ n (9), shown in refs. 0-|5[ to have the central charge of 
the g^'xg^'/g^ 2 ' conformal field theory as the limiting (y — > 00) value of c(y). Of course, 
y must be less than 9q for the above treatment to be valid, but if 9q is taken large enough 
then the relevant central charge is approached well before the approximations break down. 

At this stage it would be possible to continue by discussing what happens as y increases 
beyond 9q/2, the values of the function L a {9) near y starting to influence those at —y, and 
so on. Such a treatment would mimic that already given by Al. B. Zamolodchikov in ||18|| . 
However, at least if the aim is merely to uncover the sequence of central charges visited by 
the trajectory, it seems better to think of the above discussion just as a warm-up exercise, 
and now to proceed directly to the general situation. 

An assumption is needed, generalising already-observed behaviour of 'local' TBA equa- 
tions, on the form of L' a {9). For 9 <C —y, L and hence also L' are approximately zero, due 
to the dominance of the energy term in ( |3.5| ). Then after a transition around 9 = —y (a 
kink) during which V 7^ 0, L should settle down to some (not necessarily zero) constant 
value, and V return to zero. In the case of a local TBA equation, this is all that happens 
until 9 = +y is reached, at which point there is another kink and L returns to zero. How- 
ever, in the case of ( |3.5|) , the non-local terms in the convolution serve to link the values 
of L and hence V separated at a distance 9q, and the 'seed' kink at 9 = —y would be 
expected to propagate in and cause further kinks at — y+9o, —y+29o, and so on. Similarly, 
the kink at 9 = +y results in other kinks at 9 = y—9o, y—29o, propagating in from the 
right. This process does not form any kinks for \9\ > y since in any event the convolutions 
are swamped by the energy term in this region. 

We can try to make this a little more precise by taking the derivative of (|3.5|) with 
respect to 9, and substituting e' = — (l+e e )L' . This gives an equation which couples the 
values of V at 9, 9+9q and 9—9q. Repeating this at 9+9q, 9+29q, and so on couples L'{9) 
with L'(9+n9o) for any n G Z; furthermore so long as we remain in the region |#+n#o| < y, 
the inhomogeneous terms (2fh a e~ y smh(9+n9o) ) can be dropped. In addition, if l^+n^ol > 
y, then L'(9+n9o) = is forced (up to exponentially small corrections) by (|3.6|) . Thus 
for 9 56 ±(y—n9o) (n G Z), we expect L'(9) ~ 0. The picture is more complicated at 
the 'transition regions' around 9 = ±y, where the energy term and the convolutions are 
in equal competition. The resulting non-zero values for L'(9) at ±y then propagate in to 
— y + 9 Q , y — 9o and so forth. However, for this to be any more than a plausibilty argument, 
it would have to be shown that the set of equations coupling together the values of V at 
intervals of 9q were non-singular in a suitable sense, and also that the ^-dependent factors 

10 



of — (1 + e e ) did not spoil this. Nevertheless, it provides some support for the claim that 
L'{9) is close to zero for 9 not equal to ±(y— ndo), with n = 0,1... and y — n$o > —y, 
and that L(9) is thus a series of plateaux, separated by kinks at these points. (Note, the 
kink regions should each have a size of order one, and so taking y large while preserving 
the value of y/9o enables their size relative to y and 9q to be made as small as desired.) 
In the generic situation, kinks propagating in from the left will miss those propagating in 
from the right, the two sets of kinks being interlaced along the #-axis. This only fails to be 
true if y/9g is near an integer or half-integer; this should correspond to a crossover in the 
critical behaviour. Otherwise, take y to lie between (k— 1)9q/2 and k9o/2 for some integer 
k. Working in from 9 = — oo, the first kink to be encountered is the seed kink at 9 = —y; 
then a descendant of the rightmost kink is found at 9 = y — (k— 1)9q; then a descendant of 
the leftmost kink at —y + 0o; and so on. Between these 2k kinks, there are 2k— 1 plateaux, 
within each of which the L a (9), and hence also the e a (9), are approximately constant. The 
i th plateau, lying between the i th and (i+l) th kinks, is centred on Zi = (i—k)9o/2, and 
the approximately constant value of e a {9) along this plateau will be denoted e % a = e a (zi). 
These constants are easily found, since constant terms can be pulled out of the convolutions 
in (|3.5|), leaving just the overall integrals of </> mm and ip as proportionality factors. From 



[H and jig] , these are given by 

— / d9<P2 n (0)=S ab -2C; b 1 ; — / d9^ ab {9) = -C~ b \ (3.7) 

where C a b is the Cartan matrix of the non-affine Dynkin diagram for g. (These formulae 
can also be proved directly from ( |2.7fc| ) and ( |2.8| ) f25[|.) Considering equation ( |3.5| ) for 
9 = Zi results in the following consistency condition for the numbers e^ : 

r 

4 = J2 [{^at - M log(l+e- e *0 - C~ h x (log(l+e- £ r 2 ) + log(l+e-^ +2 ))] . (3.8) 

If Zi lies towards an end of the chain, that is if i = 1, 2, 2k— 2 or 2k— 1, the pseudo-energy 
there interacts with e^(9), already fixed by (|3.6| ); to the approximation to which we are 
working here (and increasing y while keeping y/9o fixed if necessary to sharpen things 
up), these functions can be taken to be infinite at the points at which they need to be 
evaluated. Thus, with the boundary conditions that e~ l = e° = e 2k = e 2k+1 = oo, ( |3~8l ) 
should hold for all i = 1 . . .2k— 1. Note that e l interacts with e : rather than with its 
immediate neighbours e l±1 , a consequence of (|3.5|) and the fact that ^±6*0 = ^j±2, not 
Zi±\. Thus the constants {e^} divide into two sets, {e 2 J} and {e^ J_1 }, values in each set 
only interacting with others in the same set or, at the ends, with the (effectively infinite) 
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boundary values determined by (|3.6|). This can be used to put (|3.8|) into a more familiar 
form: setting x J a = exp(— e^ 7 ' -1 ) and y J a = exp(— e^ ), the consistency conditions become 



r 



5ab—'2C r , / • -, \ C_. / ■ ,,\C. i 



6=1 



n(i+*o (i+-r) ( i+ • ■'= i - fc 

yi = IT (1 + ^f^' I 1 + d-f :h ' (1 + *r) Ci ' ■ *=*•••* 

6=1 

with the understanding that x° = x^ +1 = y\ = Ua +1 — 0. In the Ai case, there is just 
one particle type and (7-Q = 1/2; the system reduces to that given in ||, with x\ and y\ 
becoming X{ and yi in Zamolodchikov's notations. On the other hand, ( |3.9| ) is also the 
specialisation to I — 1 of the constraints for the g( k >xg( l '/g( k+l > TBA system proposed in 
|T(J: to see this, x l a , y l a above should be replaced by Y i a (— 00), Y i a _ 1 (+oo) respectively. Of 
course, in the two above-mentioned papers the 'morality' is rather different from here, in 
that (for fixed particle type a) x J a and y^ are the asymptotic values at ±00 of a number 
of different functions, indexed by j; here, they are the intermediate values of a single 
function, exp(— e a {9)). This changes the spirit of the central charge calculation, although 
the underlying structure is the same. 

One more piece of notation is needed. The interleaving of plateaux and kinks described 
above can be summarised in the following inequalities: 

z < -y < z\ < y-(k-l)9 <z 2 < -y+9 < z 3 < . . . 

(3.10) 
. . . < z 2k -i < -y+(k-l)0 o < z 2k 

(recall, Z{ = (i—k)6 /2, and (k— 1)6> /2 < y < k6 /2). Define Ki to be the interval 
[z{-%, Zi], so that the i tli kink of L a {6) lies in Ki. (For % ^ 1, 2k this is also a kink for e a , 
while for i = 1 or 2k it marks the transition of e a to the exponentially growing behaviour 
described by ( |3.6| ).) For y/9o sufficiently far from the crossover (integer or half-integer) 
values, the kinks occur well inside the regions Ki, and by the time the boundaries of these 
regions are reached, the functions L a have settled down to constants. Were it not for the 
non-local nature of ( |3.5|) , these kinks would be described by a set of independent TBA 
equations, with boundary conditions, effectively at ±00, determined by (|3.9| ). However, 
the nonlocality couples the equation for the kink in K{ to those for the kinks in Ki± 2 ; 



to leading order these coupled equations are exactly those proposed in |T(|. In this way, 
the 'tower' structure of the TBA systems in that paper can be traced back to the block 
structure of the affine Toda S-matrices, the key relation being ( |2.15| ) . 
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This would be enough to verify our claims; however for completeness the remainder 
of section concludes the computation of the function c(y, 9 ) directly in the context of the 
staircase model, continuing to work in the asymptotic region where both y and 9q are 
large, with (fc-l)/2 < y/6 < k/2. 

Equations ( [3.3|) and ( |3.5|) must be re-examined with an eye to their asymptotic be- 
haviour. The integration range [— oo, oo] in the expression (|3.3|) for the effective central 
charge c(y,6o) can be restricted to [21,22k] = Uf^-R^, since the double-exponential decay 
of L a (6) causes the integrand to be vanishingly small outside this region. Within this 

region, 

{fh a e~ y ~ 6 , 9 E K\ ; 

0/ 9eK 2 [..K 2k _ 1 ; (3.11) 

m a e-y+ e , 8eK 2k . 

Thus the only direct contribution to c(y, 9q) is from the kinks in K\ and K 2k ■ Writing 
c(y, $o) = C- + c_|_, where 

c± = —J2 d9m a e- y±e L a {e) , (3.12) 

and noting by symmetry that c_|_ = c_ , all that remains is to calculate 



TV 



V / d9m a e-y- 9 L a (9). (3.13) 

o=l Jk ± 



This is essentially standard, the only modifications needed to the cases in [^,|T5|,|T^] being, 
crudely speaking, that integration by parts should be replaced by considerations based on 
the 'locality' of </> mm (#) and i()(0). The relevant property is the following: let f(6) be any 
even function of 9 with its support concentrated at the origin (by this we mean that the 
overall integral of f(9) receives a significant contribution only in a region near #=0, the 
size of which is asymptotically small in comparison with 9q and y; both </> mm (#) and i/j(9) 
satisfy this condition). If A(9) and B(6) vary slowly (on the scale of the support of /) at 
the ends of the interval Ki = [#»-i, zi\, then 

Mf'*A(9)B(9) = - f Mf'*B(9)A(9) + ( f°° d9f(9)) [A(9)B(9)] Z ; . (3.14) 

Ki JKi \J-oo J 

If f(9) = 5(9), this formula really is just integration by parts; a proof of its more general 
validity is relegated to appendix A. 

Returning now to the calculation of the effective central charge, the necessary property 
of the solutions of the TBA equations, allowing for the exact evaluation of (|3.12| ) in the 
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asymptotic region of interest, is found by taking the derivative of (|3.5| ), multiplying by 
L a (9), integrating J K d9 and finally summing on a. In an abbreviated notation, this 
yields 



TV 2 



where 



T = irC l -T +Xl+T + (3.15) 



ji = J2 dee MLa(6) = J2 l'\de\og(l + e- £ ) ; (3.16a) 

a=l J K i a=l •'^ 

C l = 4t V / d62fh a e- y sinh(0)L a (0) ; (3.166) 

J o = ^rE / d6J>T'*L b (e)L a (d); (3.16c) 

J ± = ^E / d0fl;' ab *L b (0±e o )L a (0). (3.16d) 



These terms can be simplified in various ways. First, using the relations 

{—fh a e~ y ~ e , 9 G K\ ; 
0, deK 2 , ...K 2k _ l - (3.17) 

m a e-^+ e , 0GK 2fc , 

and referring back to ( |3.12|) , we have 

C ~ -o M c_ + <5 2M c + . (3.18) 

(Note how the relations ( |3.11|) and ( |3.17|) , by forcing the effective energy terms to lie at 
the end of the A-type Dynkin diagram, automatically select the 1=1 cases among the TBA 
systems proposed in |fT6"|.) 

Of the remaining pieces, Xq is just a surface term: applying 03.14T) to ( |3.16d ), swapping 
the sums on a and b and recalling that <^™ n/ =<^ in/ yields 

r 

2T = J2 (*«* " 2C at) [LaWWO)]*^ , (3.19) 

a,b=l 

where ([T7]) has also been used for the overall integral of </>^ m . A similar manoeuver for 
( |3.16d| ) gives, for 1 < i < 2k, 

r 

4 = ~ r T 2 - E c£[L a {e)L h {e±o Q )\*_ i , (3.20) 

a, 6=1 
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while Z^_ ~ if i < 1 or i > 2k, owing to the double-exponential decay of L a {6). The 
surface piece makes this equation a little hard to use; to get rid of it, define 



r ±= r ± + l E ^[LMUie^eo)] 2 ;^. (3.21: 

a,b=l 

Recalling that Zi+0Q = Zi+2, (|3.20|) is now 



T ± = -Xf l , (3.22) 

and, as before, T l ± ~ if i < 1 or i > 2k. Equation (|3.15[) becomes 



Z 1 

T = L-C+T_+T + + - £ [L a {e){{2C- b 1 -5 ab )L h {e) 

a, 6=1 ' 1lZ . (3.23) 

-c; b 1 (L b (e-e ) + L b (e+e )) 



Zi-l 



equation ( |3.19|) having been used to substitute for Xq. Now the term in curly brackets is 
familiar from ( |3.8| ) , since it is only ever evaluated when 9 is equal to one of the zi : once 
the sum on b has been performed, it is just e a (9). So, if we set 



3 l = T -\Y\L a {9)eMT:^ (3-24) 

and use ( |3.18|) as well, the result is 



2 ^r^v^^v^j^ 

o=l 



T = ¥ (" 5l '* c - + 6 ^ C +) +I-+I l +- (3-25) 

At i = 1, this is 

c_ = -^(J 1 -i 1 _-l 1 + ). (3.26) 

Now Jl = -J+ 1 = 0, while Z\ = -T?_. In addition, for 1 < i < 2k, ( |3T25| ) can be 
rewritten, via ( |3.22| ), as X % _ = J" 1 + X!_ +2 . Thus ( |3.26|) falls out into a cascade of terms, 
one for every other kink starting with the left-most: 



TT 2 



^{J X + J*+%L) 



n 2 



3 '' 



(3.27) 



= -^E-? 2 '- 1 



7T 
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the process terminating when the term J' 2k + 1 = is encountered. Note that the last line 
only depends on y and 9q via the value of k that their ratio entails. To calculate J 3 ~ , 
we write it as Yla=i «^a J_1 ' where 



Jl 3 ~ l = Ja 3 ' 1 -\[La{e)e a (9)] 



Z2J-1 
Z2j~2 



and 



Jl j - 



.2J-1 



' e 

= Li 



2j-2 



efelog(l + e £ ) 
■Li 2 



(3.28) 



(3.29) 



,2J-1 



r 2j-2 



Here, ~L\.2{z) = J <i£log(l — t)/t is the 'usual' dilogarithm function (for more information, 
see ref. |27] ). The surface part in (|3.28| ) turns this into the Rogers dilogarithm, C(z) = 
_ i £ dt (l2Eilzll + *£i) , V i a t he relation 



£ 



x 



1 + x 



-Li 5 



-x) 



log (1+x) log a; 



valid for x > — 1 ||271| . With the earlier definitions x 3 a 



,2j-l 



,2j-2 



( p.28|) becomes 



^ 2i - 



-C 



x{ 



+ C 



vi 



(3.30) 



equation 



(3.31) 



1 + xiJ \ 1 + yi . 
The effective central charge is c_+c_|_ = 2c_. Subject to the hypotheses on the 
plateau-like behaviour of L a (9), we have now shown that 



lim 



V,Vq— »a 

(fe-i)<2a/e <fe 



c(y,0o) 



Ck 



6 

7T 



2 EE 



£ 



V 



l + xi 



C 



ui 



i + yi 



(3.32) 



where x 3 a and y^ are the real solutions of (|3J|) (or zero in the case of y\ , corresponding to 
the first kink). In taking this limit, the ratio y/0o is assumed to be kept fixed. If instead 9q 
is fixed at some (suitably large) value and y is increased, c(y, 9q) will run through the values 
Cfc in turn. The plateau structure should break down, and a cross-over occur, whenever 
the interval between a pair of kinks is of the same order as the kink size. The intervals 
between the kinks are of lengths 2y — (k— 1)9q and k9o — 2y, while the kinks themselves 
have an extent of order one. Thus the cross-overs should occur while 2y is in a region of 
order one about each integer multiple of #o- (Examining the graphs in ref. HT8"| , the correct 
size would seem to be about five, a value completely consistent with the form ( |3.4| ) of the 
functions <j)/ x \ . The important point is that the sizes of these cross-over regions depend 
neither on 9q nor on y.) 
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To find numerical values for the c&, we use the following sum rule [28], which holds 
whenever the constants x\ satisfy a 'chain' of relations of the form (|3.9| ), of length k: 



EE 

3 = 1 0=1 



c 



X' 



l + X J a 



n- 



rhk 



6 h + k + 1 



(3.33) 



(So for the constants y l a in (|3.9| ), k should be replaced by k— 1 here.) For the simply-laced 
Lie algebras, the Coxeter number h is equal to the dual Coxeter number h, so the right- 

2|] at 



hand side of ( |3.33| ) is in fact the central charge for a theory of Gepner parafermions 
level k+1. Since we do not know of any general proof of these relations, we checked them 
numerically (to eight-digit final precision) for all simply-laced Lie algebras of rank < 8, for 
k — 1 . . . 8. With the relation 



C 



x 



1 + x 



~6~ 



this gives 



c fc = r(h + 1) 



1 



C 



k 



1 



1 + x 



fc + 1 



h + 1 h + k h + k + 1 
a result which reproduces the central charge of the g( k >xg( 1 >/g( k+1 > coset model. 



4. Conclusions and speculations 

The hopping nature of the solutions to equation ( |3.5| ) has now been established in a fair 
amount of detail, but the physical meaning of this result is much less clear. While we have 
no final answers to this question, we give in this concluding section some speculations on 
the Lagrangian form of the staircase models which might be worthy of further investigation. 

The usual Lagrangian density for affine Toda field theory can be written as 



C 



1 



(dcf>) 2 -V(^(3) 



where 



m 

J 2 



n'V'^^EV**-*. 



(4.1a) 



(4-16) 



i=0 



the sum in the second equation running over all the simple roots a.\ . . . a r of (non-affine) 
g, together with a — — Xa 71 * ^ the negative of the highest root. So long as (3 is real, 
the potential has a unique minimum at (p = 0, and a purely elastic scattering theory of r 



different types of bosonic particles results [p0| -p3 . 
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However the potential ( }4.1£| ) has also been discussed at purely imaginary values of (3 
, the main aim being to find a Lagrangian basis for the 4>ud,id,adj) -perturbed coset 
models (generally thinking of the perturbation as being in the massive direction). In such 
cases the unitarity of ( |4.1a|) is doubtful, and some form of truncation has to be invoked to 
restrict the space of states to a positive-definite subspace. The physical motivation for this 
procedure is unclear, and so it might be preferable to look for an alternative continuation 
of the real coupling-constant Lagrangian for which the manifest reality of the potential 
was preserved. 

The simplest way to create a real potential from ( f4.1fc|) is just to add its complex 
conjugate, forming V '((/), (3) = V(<j>,/3) + V(4>,(3*). (We assume that the field <fi remains 
real here.) However, the non-linearity of the To da field equations means that at the 
classical level there is no reason to believe the potential V to be integrable. But at the 
quantum level, for suitably well-chosen values of f3, things may be different. The basis 
of this expectation is the strong-weak coupling duality of the affine Toda field theories, a 
quantum effect seen for the S-matrices in the fact that (from equation ( |2.2|) ) S^ b = S ab . 
Since B(/3) — 2(3 2 /{(3 2 + 47r), the interchange between the two possibilities is effected by 
the following 'duality' transformation on the coupling constant: 



(3 - = 471-//?, (4.2) 

relating the strong and weak coupling regimes. At the more fundamental level of light- 
cone canonical quantization, it has also been argued p3] that the Hamiltonians arising 



from V((j),P) and V{(j),f3) actually commute - a result from which the equality of the 
strong- and weak- coupling S-matrices follows automatically. This motivates the idea that 
a linear combination of V(<p, (3) and V((p, (3) might be quantum integrable. One possibility 
is to form the 'self-dual' potential V ,f (<f>,/3) = V((f), 0) + V(4>,(3), and in fact analogous 
non-affine potentials have already been examined by Mansfield [|35|| , though with rather 



different motivations. Here, the hope is to construct a real continuation of Q4.1|) which 
retains integrability at the quantum level, and we now see that this might be possible if 
V'((/), [3) = V"((j), 13), that is if (3* = 0. The allowed values of (3 are thus 

j3 = 2^e l \ (4.3) 

and the 'roaming' analytic continuation has been recovered, from an alternative perspec- 
tive. As a function of the (real) variable 7, the proposed potential is 

2 r 
V"(</>, 7 ) = — V nie 2v^cos 7a ^ cog r(2 v ^Fsin 7 )a i .(/) + 2 7 1 . (4.4) 

27T L — ' 
i=0 
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Taking 7 = recovers the real-coupling Toda theory at its self-dual point, while (since 
9q = j-taxi'y) the large-#o limit corresponds to 7 — ► n/2. In this limit each term in the 
sum ( [4.4] ) is a product of a slowly-growing exponential with a rapidly-varying cosine, and 
near the origin the total is well approximated by 

2 r 
--— y^nicos r^v 7 ^ sin 7) «;.</> + 27] . (4.5) 

Z7T z — ' 
i=0 

We would like to think of this as the 'unrestricted' potential, a manifestly real analogue, 
for the higher-rank algebras, of the sine-Gordon potential. In the A2 case, for example, 
the potential ( f4.5| ) has minima lying at the vertices of a triangular lattice. If these are 
the classical vacua, then (at least from a visual examination of the potential) kinks should 
exist, joining neighbouring vacua, one for each edge of the lattice. This motivates the idea 



that (|4.5|) might provide a potential for a Landau- Ginzburg description of an unrestricted 
SOS-type model associated with the algebra g. Further from the origin, the exponential 
growth of the other terms in ( |4 . 4| ) becomes important. At first sight this is unfortunate, 



as an immediate consequence is that the potential ( |4.4| ) is not bounded from below, and 
so it seems that the theory is not well-defined. However (at least so long as #0 is large 
enough) while the minima get lower at larger values of |0|, so the saddle-points between 
these minima get higher. At a given energy-scale, this may serve to confine the field to 
a neighbourhood of the origin, the vacua further away being effectively decoupled. To 
take a classical analogy, a particle starting at the origin of the potential ( |4.4| ) with any 
finite kinetic energy would be confined to a finite domain in the emplane, despite the fact 
that beyond this domain points can be found at which the potential energy is arbitrarily 
negative. The question is thus whether 1+1 dimensions is near enough to the classical 
(d = 00) limit. Assuming that this is true, a rather appealing picture emerges - - at 
any given energy-scale, only a certain number of vacua are accessible to the field (j), and 



the SOS picture provided by (|4.5|) is modified to an effective Landau-Ginzburg potential 



appropriate for a restricted SOS-type model (the possibility of such potentials for the case 



of ct2 is mentioned in [10] ). As one progresses towards higher energies in the UV limit, 
more and more vacua become available and the restriction is lifted in stages, in a way which 
could reproduce the roaming behaviour already observed at the level of TBA equations. 

It has to be admitted that there are many problems with this scenario, even at the 
naive classical level presented above. In particular, it seems hard to ensure that the easing 
of the restriction happens in the correct manner, so that a full series of coset models 
(or rather their hoped-for Landau-Ginzburg potentials) is visited. In this regard, it may 
be relevant that there is no compelling reason for the choice of V" given above - more 
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generally, pick any function /, with complex conjugate /. Then f((3)V((f), (3) + f((3)V((f), (5) 
is a candidate potential. Clearly it would be more satisfying to have a physical principle 
to justify one choice over all the others. It should also be stressed again that the models 
described by (fOj) are perforce a long way from their classical domains, and so (just as 
for the Landau- Ginzburg picture of the c < 1 models |36| ) classical intuitions should be 
treated with caution. 

Despite these difficulties, the search for a physical understanding of the staircase 
models seems to be very worthwhile. At the level of the TBA, they already unify and 
clarify a whole class of systems; it is not unreasonable to hope for similar insights at the 
field theory level as well. 
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Appendix A. 

The aim is to show that 

d0f'*A(O)B(O) = - [ <Wf'*B(0)A(0) + ( r Mf(6)) [A(0)B(0)]* , O 

'Ki JKi \J-oo / 

subject to the conditions already mentioned in section 3. First, expand out the left-hand 
side as follows: 

d9f'*A(0)B(0)= f del de'f'{9-e')A(0')B{9) 

Ki J Ki J — oo ( A ~\) 

d0de'f{d-6')A(0')B{6), 

IUIIIUV 

where I, III and V are three regions in the 9, 6' plane: 

I = Ki x [-oo, zj_i] ; III = KiX [z i: oo] ; V = Ki x K % . 



Similarly, the first integral on the right-hand side of ( ^.14| ) is, with a trivial change of 
integration variable, 

de'f'*B($')A(e') = - I d$' I d9f'{9' - 6)B(6)A(6') 

Ki JKi J — oo /a y\ 

dede'f{e-6')A{e')B(e), 

nurvuv 
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swapping the order of the integrations (the double integral is absolutely convergent due to 
the finite extent of Ki) and using the antisymmetry of /' in the last line. This time the 
last integral is over the three regions II, IV and V where V is as before and 

II = [-00, Zi-i] x Ki; IV = [zi, 00] x Ki . 



Subtracting ( |A.2| ) from ( |A.1| ), the second term on the right-hand side of ( |3.14j ) must be 
given by 

d9d9'f(9-9')A(9')B{9). 



+ 
L<T <TI "'III ./IV J 

Consider first the integrals over I and II. The locality of /' means that in these two regions 
the integrand is only non-negligible near (zi—i,Zi-i), and the approximately constant na- 
ture of A(9')B(9) near this point allows its replacement in the integrand by A(zi-i)B(zi-i). 
Asymptotically we are also free to replace I by I' = [zi-i, 00] x[— 00, z$_i] and II by 
II' = [— 00, Zi—i]x[z{-i, 00] (it may help to draw a picture of this), after which the re- 
maining (double) integral of f{6 — 9') is easily seen to reduce to — J f(9). In just the 
same way the difference between the integrals over III and IV is well-approximated by 
(f_ f(6))A(zi)B(zi), which completes the demonstration of ( |3.14| ). 

Note that in the usual TBA, the integrals all run from —00 to 00, in which case a 
naive swap of integration variables leads to the second term on the right-hand side of 
(J3.14J) being missed. This would be wrong, the reason being that the integral over the 
entire 9, #'-plane is not absolutely convergent. The correct prescription is to regulate the 
6>-integral to an interval Ki, say, and then let Ki — > [—00, 00] at the end, thus recovering 
( P-14J ) with limits of ±00. It is worth stressing that this result does not depend at all on 
the 'integrable' nature of the TBA equations, equations which only need to be used once 
(in the form of the relation ( |3.15|) ) at each step in the evaluation of the effective central 
charge. 
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